Ideal and nonideal electromagnetic cloaks 
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We employ the analytical results for the spatial transformation of the electromagnetic fields to 
obtain and analyze explicit expressions for the structure of the electromagnetic fields in invisibility 
cloaks, beam splitters, and field concentrators. We study the efficiency of nonideal electromagnetic 
cloaks and discuss the effect of scattering losses on the cloak invisibility. 



Cloaking of objects was first explicitly suggested by 
Pendry Q as a method to cover an object by a layer of 
a composite material with varying characteristics in such 
a way that the scattering of electromagnetic waves from 
this object vanishes, so that the structure of the electro- 
magnetic field outside the cloak appears as if the object 
and cloak are absent. The required material for creat- 
ing such invisibility cloaks should have rather complex 
anisotropic, spatially varying properties. Moreover, such 
a material should possess both nontrivial dielectric and 
magnetic responses. Generally speaking, it is impossible 
to find dielectric media with the required varying prop- 
erties in nature, and therefore the first electromagnetic 
cloak realized experimentally employed microwave meta- 
materials Q. Metamaterials are artificial microstruc- 
tured materials with electric and magnetic properties de- 
signed by careful engineering of their constituents. The 
constituents of these composite materials are normally 
sub-wavelength electric and magnetic resonators, with 
the specific properties determined by their geometry. The 
metamaterials can also be anisotropic, and thus they can 
satisfy, in principle, all the requirements for realizing in- 
visibility cloaks for electromagnetic waves. 

The concept of the electromagnetic cloaking is based 
on the coordinate transformation of space and the cor- 
responding transformation of Maxwell's equations, which 
provide the required expressions for the effective spatially 
varying dielectric permittivity and magnetic permeabil- 
ity of the cloak medium The possibility of creating 
such cloaks has been confirmed experimentally and 
it was further discussed and demonstrated in many nu- 
merical studies [3, i, i, i, [i, ,8,, 1, M] ■ Numerical simu- 
lations are usually performed either by using commercial 
finite-element equation solvers P, 0, H, H, or by em- 
ploying the decomposition of the electromagnetic fields 
into a set of the corresponding eigcnmodes [1, 0, S @1 ■ 
Several studies have analyzed a possibility of creating 
simplified cloaks 0, , since the required media parame- 
ters for an ideal cloaking are extremely complex for their 
realization in experiment. In particular, an optimiza- 
tion of the cloaking parameters for suppressing scattering 
has been analyzed in Refs. 0, 0- In this context, it is 
worth mentioning an earlier pioneering study [ll[ , where 
an analogy of the coordinate transformation to an effec- 
tive magneto-dielectric medium was discussed for some 



particular cases. 

In this Letter we analyze the explicit analytical ex- 
pressions for the structure of the effective electromag- 
netic fields caused by the spatially varying coordinate 
transformations. In a sharp contrast to the rather formal 
transformations of Maxwell's equations in the Minkovskii 
form used in the original pioneering paper [l| and 
all subsequent publications, we present a simpler form 
of the time-independent transformations that can be im- 
mediately employed for solving many important prob- 
lems. Based on these results we suggest a simple recipe 
for evaluating the screening efficiency of non-ideal cloaks 
and calculate their scattering dipole moment. 

General formalism. We start our analysis by consider- 
ing the propagation of monochromatic (i.e. ~ exp(iwf)) 
electromagnetic fields described by Maxwell's equations 



curl H 
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-B 



(1) 



First, Eq. ([T]) can be written in the component notations 
in the following form. 



ox" c 



(2) 



where E is the totally antisymmetric symbol [13[, 
is the orthonormal basis corresponding to the Carte- 
sian coordinates, the indices (a, b, c) run through the val- 
ues (1,2,3), and are the Cartesian coordinates corre- 
sponding to the orthonormal basis . In a medium with 
the dielectric permittivity ig, the electric displacement 
vector D can be expressed in terms of the electric field 
as = i^^Eb. The magnetic displacement vector B can 
be expressed as i?" = in a medium with magnetic 

permeability fiQ. We note that the repeated upper and 
lower indices imply summation (Einstein convention). 

Now we introduce new coordinates y through the 
transformations y = y(x) with the metric tensor g^^ de- 
fined in a standard way, 



dx'^ dx^ 



(3) 



where 5°'^ is Kronecker's delta-function. The determinant 
g{x) = det ||5*-'(x)|| is of a particular importance for the 
coordinate transformations. 
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In the framework of these new coordinates the Lh.s. of 
Eq. ^ can be re-written in the form, 



2 \dx'' dy^ dx''- dy^ 



l^abc I dy^dy^dHk _ dy^ dy' dH, 
2 I dx^ dx'^ dy^ dx'^ dx^ dy^ 



(4) 



where H is the transformed magnetic field defined as 



~ Ox^ 
oy'^ 



(5) 



The totaUy antisymmetric symbol is a tensor density, 
and real tensor (Levi-Civita tensor, e*^*^), which coincides 
with E'^^^ in Cartesian coordinates, is related to the E"^^^ 
through e^^^ = ^E^^'^, and it transforms under coordi- 
nate transformation according to the tensor rules 



Jik 
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dx'^ dx'^ dx'^ 
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Multiplying both parts of Eq. ^ by dy'' /dx"^ and using 
Eq. ([6]), we obtain 
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where E is the transformed electric field defined as 

En - g^E, , E, = —^E^ . (7) 
As a result, we arrive at the transformed equation 



E - (cuilH) - —^T^T^Sq En , (8) 



dy^ 



c dx°- dx^ 



which coincides with the original Maxwell's equations but 
it is written for the new electromagnetic fields in the new 
(spatially transformed) space. The left-hand side of this 
equation contains the curl operation on the transformed 
magnetic field, whereas in the right-hand side of Eq. Q 
we have a new electric displacement vector which is de- 
fined as the product of the new dielectric tensor, 



J_dT/_djf_ „j 
y/gdx'' dxb^° 



(9) 



and the new electric field ([7]). Applying the same trans- 
formation procedure to the second equation of Eq. ((1]), 
we obtain the corresponding expression for the effective 
magnetic permeability. 



,„ _ 1 dy^ ay" 



^ dx°- dx'' 



(10) 



When the initial medium represents vacuum (eg" = 
^g" = J""), the resulting new dielectric tensor coincides 
with the effective magnetic tensor. 



Our results for the effective material parameters coin- 
cide with those of Ref. Q but, in addition, they also 
provide explicit analytical formulas for the transformed 
fields. Thus, we can find analytically the field distribu- 
tion for any specific cloak without solving Maxwell's equa- 
tions in the transformed geometry, but just transforming 
the free-space fields. We have both the expressions for 
the structure of the cloak and equations (O and ([7]) for 
the field transformations in the cloak. As a straightfor- 
ward application, our results allow analytical estimations 
of the scattering from non-ideal cloaks, and thus a po- 
tential optimization of their parameters. 

Examples. We use those results for calculating the 
field structure in both cylindrical and spherical cloaks, 
field concentrators, and other devices designed via the 
transformation optics. As the first example, we consider 
the linear coordinate transformation for the cylindrical 
cloak [11, X = X{x, y),Y = Y{x, y), Z ^ z, where the 
radius is transformed as i? = a + r{h — a)/h, where R = 
vA^ + y^, r ~ ^ x^ +2/^- The cloak occupies the space 
a < R < b, where a and b are the inner and outer radii of 
the cloak, respectively. To find the field distribution in 
the cloak illuminated by a plane electromagnetic wave, 
we write the plane wave in vacuum as = exp(i«;yy), 
Hy ~ 0, Ez = cxp{iKyy) . Applying now the transforma- 
tions according to Eqs. ([5]), we find the field distribution 
in the electromagnetic cloak with the material parame- 
ters © and (Uni), Ez ^ cxp{iKyY)[b {b - a)~^ {1-aR-^)], 
= b{b - a)-i(l - aY^R-^)Ez; Hy = abXY{b - 
a)^^R^^Ez, for a < R < b. The field vanishes inside 
the cloaked area, whereas outside of the cloak the field 
is an unperturbed plane wave. At the external surface of 
the cloak (at R = b), the normal component of the mag- 
netic field is discontinuous, since the linear space trans- 
formation produces discontinuities of the cloak material 
parameters. Using the transformation function 



6- i? 



(11) 



provides, for large enough values of /3 (/? > 1) continu- 
ity of the field components at the external surface of the 
cloak. As a matter of fact, Eqs. (O and ([7]) allow us to 
find the field structure for an arbitrary continuous space 
transformation. Though for an arbitrary transformation 
one cannot find the analytical expressions for the elec- 
tromagnetic fields, simple numerical calculations provide 
the required results without solving Maxwell's equations. 

As an example of a cloak of complex shape, we con- 
sider cut-type space transformation where the cloak cov- 
ers an arbitrary area 'cut' from the space (see, e.g.. 
Fig. 1). We define the internal surface of the cloak 
as a set of all points which are closer than b to the 
cut line Xcut(C)) ^ut(C)- For each point {X,Y) we 
find the corresponding closest point on the cut line 
{X^^^,Y^^^.)(X,Y) by finding the minima of the func- 
tion [(A - Aeut(e))' + {Y- Fcut(O)']'/' for all possible 
^. We define the distance from the 'cut' as R{X,Y) = 
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FIG. 1: (color online) Structure of the magnetic field Hx for 
(a) an arbitrary-shaped cloak and (b) multiple cloaks pro- 
duced by the space transformation in accord with Eq. 
where we take the parameters a — 0.3, b = 0.6, and 13 — 2. 
The cloak is illuminated by a plane wave, H° = exp(iKj,j/) for 
Hy = 5.4. Center of the cloak is shown by red line, in the case 
(a), and dots, in the case (b). 



[{X - + (Y - y'e"ut)2]^/2 and create the cloaking 

area a < R < b hy transforming the space according to 
Eq. pT|) . In the area of the cloak, a < R < b, we define 
the new coordinates [x, y) as 

X = + {X- XO„J^ , y = Y^^, + {Y- , 



R 



R 



which provide a conformal mapping of (x, y) into (X, Y). 

Figure 1(a) shows the distribution of the x-component 
of the magnetic field for the plane wave incident on the 
cloak. The cut is created through the graphic input in 
Matlab, where we specify several points (five points, in 
this example), which are connected by a smooth line 
through the cubic spline function. Figure 1 (b) shows the 
field distribution in the cloaking problem with several 
cylindrical cloaks centered at the same five points. 

Parameters of transformation r{R) are a = 0.3, b = 
0.6, /3 = 2. The field distribution is calculated for the 
incident plane wave Hx = exp{iKyy) with Ky = 5.4. 

We note that, despite the divergence of the dielectric 
and magnetic functions at the internal interface of the 
cloak, the electric and magnetic fields remain finite. For 
example, the radial and tangential components of the 
magnetic field can be found as 



\X\ 



(b-a) R 



\Hr\ = 



\Y\ 



{b-a) R 



We note that \Hfi \ is finite at the internal interface of the 
cloak, while \Ht\ vanishes when R ^ a. 

From the uniqueness theorem of electromagnetism we 
can conclude that since the tangential components of the 



electric or magnetic field at the closed surface vanish, the 
field in the volume surrounded by this surface vanishes as 
well; this means that we indeed have the volume {R < a) 
concealed from the electromagnetic radiation. 

Effectiveness of non-ideal cloaks. Above we demon- 
strated that a cloak with the parameters defined by 
Eqs. ([9]) and (fTO|) provides complete invisibility of the 
objects hidden inside. Moreover, such a cloak works for 
an arbitrary transformation function y = y(x), either 
smooth or piece- wise continuous [3]. However, this kind 
of ideal cloaking seems to be impossible in practice, since 
it requires infinite values of material parameters at the 
internal surface. The simplest way to overcome this dif- 
ficulty is to truncate the transformed area. For example, 
for the coordinate transformation with the internal ra- 
dius a we create the cloak only for R > a + Sa, and then 
place a perfectly conducting metal surface at i? = a + Sa. 
Now, if we make an inverse transformation from the cloak 
to vacuum, the metallic surface transforms into a metal- 
lic cylinder with the radius roff. If dR/dr\r=a 7^ 0, then 
Toff w 5a(dRl dr\r=a)~^ 1 the simplified cloak scatters light 
as a metallic cylinder of the radius res, which can be 
much smaller than the size of the hidden region a. 

One of the main limitations of the resonant medium 
used for creating electromagnetic cloaks is losses. To es- 
timate the effect of small losses on the cloak performance, 
we assume that e = jl = e' + ii", e" <^ e', whereas e' is 
given by Eq. ([9]). We consider a plane wave incident on 
such a cloak. After applying the inverse coordinate trans- 
formation, which transforms the electromagnetic fields 
into vacuum plane waves, we obtain in original space the 
excited effective electric and magnetic currents depen- 
dent on vacuum electric and magnetic fields. 



Je{r. 



ah 



jj^/g dx^ dx^ 
4tt dy^ dy^ 



These currents lead to the scattering radiation which can 
be easily evaluated either analytically or numerically. 

Next, we can estimate the effect of mismatched dielec- 
tric permittivity and magnetic permeability on the cloak 
performance. We assume that e = jl + 5e = Eq -\- 5e, 
where Eq corresponds to an ideal cloak and \&e\ <C |eo|- 
Then, we apply the inverse transformation x ^ y and 
find that in the original space the dielectric permittivity 
differs from the unity by the value 



8e 



xab 



oy^ ay^ 



If we assume that (e — /t) oc £, namely 5e^^ = asg (for 
a <C 1), then Si""^ = aS"^^ . Thus, in this case the non- 
ideal cloak scatters the waves as an object (cylinder of the 
radius b for the cylindrical cloak) with a scalar dielectric 
permittivity \ + a Ki 1 . 

Now we consider Se"^^ = aJ*^ , for a <C 1 . This gives us 
a relatively simple result for 5e°'^ ^ 



6s 



xab 



dx^ 

Qyl Qyt 



[12) 
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FIG. 2: (Color online) Scattering patterns for non-ideal cloak 
for TE (a) and TM (b) waves. Red lines (1) indicate scatter- 
ing from bare metallic cylinder, while blue line (2) corresponds 
to the scattering from the cloaked cylinder with truncation 
level ro = 0.056 and a = 0. Black curves (3) correspond to 
the scattering (mainly forward) on the dielectric cloak with 
parameter a = 0.1 for the figure (a) and a = 0.02 for the 
figure (b). 

and also an expression for the electric current 

Jeff « ^«£""'E", (13) 

which can be used for calculating the scattered fields. It 
is worth mentioning that in this case the cloak is sen- 
sitive to small perturbations of the parameters. Indeed, 
the inverse dielectric permittivity has a singularity at the 
inner surface of the cloak, so if the truncation value Sa is 
small then the parameter a has to be very small as well 
so that it does not perturb the structure of the plane 
wave. This result has been predicted earlier in Ref. [1]. 
We also note that very similar expressions for the electric 
and magnetic currents in the lossy case may also lead (for 
particular forms of e") to such a strong sensitivity. 

As an example let us consider plane wave scattering on 
a slightly non-ideal cloak with the magnetic permeability 
tensor components firr = [R ^ o,)/R^ fJ-(t>(i> ~ R/{R — a), 



jjLzz = 7~^(i? — a)/R, J = {b — a) /b and dielectric per- 
mittivity tensor Su = i^a + a, a ^ 1. We suppose 
that the wave propagates in y-direction, its wavenum- 
ber, ky = 5.4, and inner and outer radii of the cloak are 
a = 0.3, b = 0.6, respectively. For the case of a = such 
cloak is ideal, and it corresponds to the linear coordinate 
transformation R = a + 7r, < r < b. In order to avoid 
the singularity we place a metallic cylinder at _R = a + Sa, 
which conceals the inner region of the cloak but it also 
distorts scattering performance. Scattered fields are re- 
garded as weak, and we use the perturbation approach 
to calculate them. Let us make an inverse coordinate 
transformation i? — > r, so that the metal screen radius 
transforms to tq = 5a/ j. We calculate the effective cur- 
rent according to the Eq. (|13p and find the scattering field 
pattern for both TE (E" = z^E) and TM (H" z"i7) 
waves (see Figures 2a, b). One can sec that for the TM 
polarization the scattering from non-ideal cloak is negli- 
gible, and in the shown scale it is represented by a point 
in the center of the Fig. ©(b). If the truncation level, 
To, is reduced, the scattered radiation for TM wave is 
also reduced. On the contrary, TE wave scattering is 
significantly enhanced, e.g., if we decrease ro to 0.0056, 
than the scattered power becomes four times larger. 

This approach allows to estimate the scattering losses 
in an arbitrary nonideal cloak (e.g., such as that shown 
in Fig. 1), since the effective currents can be calculated 
using Eq. (jl3p with the plane wave excitation. 

In conclusion, we have employed the analytical results 
for the distribution of electromagnetic fields in the trans- 
formed coordinates to analyze an explicit solution for the 
problem of the wave scattering in electromagnetic cloaks, 
beam splitters, energy concentrators, and other devices 
of the transformation optics. Employing those results, 
we have derived a simple criterion of the efficiency of 
nonideal cloaks based on estimations of wave scattering. 
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